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Abstract 

A simple closed form expression is obtained for the scattering phase shift per- 
turbatively to any given order in effective one-dimensional problems. The re- 
sult is a hierarchical scheme, expressible in quadratures, requiring only knowl- 
edge of the zeroth order solution and the perturbation potential. 
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Since most quantum mechanical problems cannot be solved exactly, perturbation theory 
is a useful tool. Standard perturbation theory relies on the use of the Green function or 
the spectral summation over the intermediate states. Variants of the standard perturbation 
theory were introduced by Schrodinger Podolsky ||, Sternheimer || and Dalgarno and 
Lewis [|| . Originally, the Dalgarno-Lewis method was introduced as a means of calculating 
sum rules. Another useful variation of perturbation theory was introduced by Wentzel || 
and later by Price |J and by Polikanov f?j. However, this method was overlooked until its 
rediscovery by Aharonov and Au || and more or less about the same time by Turbiner ||. 
This method has now come to be widely known as "Logarithmic Perturbation Theory" or 
LPT for short. 

LPT has many advantages. In one- dimensional problems, it becomes possible to obtain 
perturbative solutions to all order in a hierarchical scheme, for the energies and wave func- 
tions for the bound states f?|-f| and for the phase shifts in scattering states flOU - A few 



modifications of LPT to handle the presence of zeros in the wave functions in excited bound 



states were given by Au et. al. JTT 



In three-dimensions, the zeros of excited bound states appear as nodal surfaces. A 
variation of LPT to address this difficulty was given by Au JT2]. In this method, the bound 
state wave function is written as Fexp(-G), and a perturbation expansion is carried out 
only on F. In the absence of any nodes, the zeroth order function F can be set equal to 
unity as the information on the unperturbed wave function can be totally absorbed into G. 
This is then equivalent to writing the perturbed wave function as a scalar function times the 
unperturbed wave function. When confined only to the first order correction, this method 
yields the equivalent of the Dalgarno-Lewis method. Henceforth, we refer this method 
of writing the perturbed wave function as a product of the unperturbed wave function 
and a scalar function as the Dalgarno-Lewis Perturbation Theory (DLPT for short). The 
difficulties associated with the presence of zeros in the wave function when using DLPT are 
also discussed in reference 1T2] . 



In Eq. (1.25) of Ref. [0], the hierarchical structure of DLPT is displayed explicitly 



for a multidimensional system. For one- dimensional system, this hierarchy can be trivially 
integrated twice to give the perturbative solutions in quadrature to reproduce the results of 
Kim and Sukhatme [pjfl . Recently, Nandi et. al. []14| applied DLPT to the scattering phase 



shifts in one-dimension and rederived the results of Au et. al. [|10j obtained via LPT. 

The equivalence of DLPT and LPT was first discussed by Au and Aharonov |I5[], and 
subsequently by many other authors |]16| . In the absence of nodal difficulties, this equivalence 



can be expected simply from the well defined mapping between the wave function and its 
logarithm. ^From the Wronskian conditions, one can easily see that the scattering wave 
function in one-dimension with the expected asymptotic behavior is free from zeros and 
hence its logarithm is regular. It follows from the above that LPT and DLPT should 
produce the same perturbative scattering phase shifts. 

In this letter, we apply DLPT to the scattering problem in one-dimension and derive a 
closed form expression for the perturbative correction to the phase shift to any order. This 
is an improvement over earlier works where the calculational procedures were laid out, but 
no such closed form results were available. 

The unperturbed wave function ipo( x ) is given by the Schroedinger equation 

- + V(x)i/> (x) = |ArVo(z), (1) 

for a particle with energy |fc 2 . In our notation h — m — 1. Similarly, ip(x), the perturbed 
wave function, is given by 

- + V(x)ip(x) + XU(x)ip(x) = \k 2 ij){x), (2) 

with XU(x) the perturbing potential. Both Schrodinger equations are subjected to the 
boundary condition at infinity. 

ipo(x) = ip(x) = exp(-ikx), as x — >• oo. (3) 

Then the phase shifts are defined as their respective phases at x — 0. 

exp(-2z5 ) = V *(0)M)(0), exp(-2iS) = ^*(0)/^(0). (4) 



The Wronskian condition, 

ip (x)ipQ (x) - $J(a;)V>o(aO = 2ik, (5) 
prevents ipo{x) from having any node, and ensures that the ratio 

f(x)=^(x)/Mx) (6) 
is well-defined. The perturbative correction, in terms of f(x), is given by 

Mx)f"{x) + W {x)f'(x) - 2XU{x)%{x)f{x) = 0. (7) 

where one can perform a perturbative expansion on f(x), 

f(x) = f (x) + Xf 1 (x) + X 2 f 2 (x) + ... : (8) 

with fo(x) = 1 identically. Performing this expansion on Eq. (|7|), one gets to zeroth order, 

A : ^(x)fZ(x) + 2%(x)f' (x)=Q (9a) 
and to the nth order, 

A" : Mx)fZ(x) + 2^{x)f' n {x)= 2U(x)Mx)fn-i(x), n>l. (9b) 

The first one vanishes identically as fo(x) = 1, while the second one can be recasted in the 
form 

^l{x)f' n {x))' = 2U{x)^ 2 (x)f n ^{x), (10) 
which can be integrated from infinity to give in a hierarchical scheme: 

roc j roc 

/„(*)=/ dy-— dz2U(z)i> 2 (z)f n -i(z). (11) 

As a result, one can apply it recursively to obtain f n (x) of arbitrarily high order. 

Relation ( |TTj) can be solved in terms of the zeroth order information. The Wronskian 
condition (|5]) can be rewritten as 

^)(tt\) =P(x)q'(x) = 2ik, (12) 

\Mx) 



where the definitions of p(x) and q(x) are identical to that in Ref. JlQ|JTT[ . 



^=W)~m- (13) 



Then relation ((TTJ) becomes, upon using these two relations: 
f n {x) = / dz2U{z) p(z) f n -i(z) / 



1 <° :x 

ik 



p{y) 

dz2U(z) p(z) /„_!(*) - = J[/„-i](x), (14) 



with 



1 



J[g}(x) = -I dz2U(z)p(z)(q(z)-q(x))g(z). (15) 



So one can apply Eq. (0) iteratively and gets 

n times 

f n (x)=7[Jl~?[l]---]](x). (16) 

After obtaining the expressions for f n (x), the only remaining task is to express 5 n in 
terms of these / n 's. Recall that the phase shift corrections 5 n are naturally defined by 

5 = 5 + \5 1 + X 2 5 2 + ■■■. (17) 

^From the definitions of f(x) and 5, it is easy to see that 

5-5 = Imlog/(0). (18) 

Performing an expansion in powers of A, one can express 5 n as a function of the n-th 
derivative of log/(0), 

^ to M^ iog/(o) L' (i9) 

Now we have obtained the correction of f(x) to arbitrary order by Eq. (|T6"D, and expressed 
S n in terms of /(0) by Eq. (|n3D. One only need to combine these two pieces of knowledge to 
get a closed form expression of S n . This can be most conveniently achieved by the lemma 
below. 



Lemma: If g = g(f(X)) and /(A) = 1 + A/i + A 2 /2 + • • •, then 



1 d n £ 



= £ ■ } 9®f£&-f!r, ( 2 o) 

A=0 {i p } n 'l'*2- " ' 



n! d\ r 

where {i p } n = {(ii,i2, ■ ■ ■ ,i n )} is the set of n non-negative integers satisfying 

n, (21) 



n 


P =i 








d/i 


/=i' 



j is defined by j = I]™ =1 i p , and g^ 

The proof is straightforward, just perform the double Taylor expansion and collect terms 
of the same order in A. For our purpose, g = log / and hence gV> = (— — 1)!. This 
leads us to the central result of this paper, 

5 n = Im £ (-l)J- 1 , \ J ~ 1)1 .../;«, (22) 

where /„ = / n (0) is given in Eq. (|TJ). 

Physical results should be independent of which perturbative scheme one chooses to 
work with, and the S n obtained above should agree with the standard Rayleigh-Schrodinger 
results as well as those from LPT. The agreement of LPT with the Rayleigh-Schrodinger 
theory on phase shifts was demonstrated explicitly in Ref. [10| up to order A 4 . We shall 
demonstrate below that Eq. (p2|) also reproduces exactly the same expressions for S n given 
by LPT. 

For any given n, we are going to enumerate all the non-negative integer n-plet 
) which satisfy the condition J2p=iPip — n - Each of these n-plet is going to 
give a term in Eq. (|22|). The integral expressions for / n (0)'s are substituted in, and the rest 
are just simplifications. 

For n = 1, there is only one element in {i p }, namely (1). Then Eq. (^) gives 

1 f°° 

Si = Im fx = Im — / dx U(x) p{x) q{x) 

%K J 
— 1 f°° 

= — — Re / dx U (x) p(x) q(x) . (23) 
k Jo 
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Note that g(0) = 0. This agrees with the LPT result (Eq. (26) in Ref. |L(J). In the notation 
defined in Ref. [HJ, with 

/*oo /*oo /*oo 

I[F u F 2 ,...,F n ] = / dx 1 F 1 (x 1 ) dx 2 F 2 (x 2 ) ... dx n F n (x n ), (24) 
the result is 

5 x = ^-ReI[Upq\. (25) 

This notation will facilitate the presentation of higher order results. 

Less trivial is the case n = 2, where (ii,i 2 ) can be either (0, 1) or (2,0). Eq. (21) then 
gives 

5 2 = Im(/ 2 -i/ 1 2 ) . (26 ) 

In particular, 



oc 



h=-nr\ dxU{x)p{x)q{x) I dy V \y) p(y) (q(y) - q(x)) 



k 2 Jo 

ll[U Pq 2 ,Up] + ^I[Upq,Upq] 

^I[Upq 2 ,Up] + lfl (27) 



As a result, we have reproduced Eq. (31) in Ref. [|10|] : 

5 2 = ^lmI[Upq 2 ,Up}. (28) 

The n = 3 is more cumbersome but still straightforward. In this case (ii,i 2 ,i3) can be 
(0, 0, 1), (1, 1, 0) and (3, 0, 0). Hence 

( 5 3 = Im(/3-A/ 2 + |A 3 ). (29) 

After expanding out the / n 's and some changes of variables, the result can be casted into 
the form 

S 3 = ^ReI[U Pq 2 ,Upq,Up}, (30) 



which has been reported in Eq. (33) of Ref. . 



Similar exercises can be done with even higher n. We do not think, however, that it 
is instructive to present those higher order calculations in full details. We just want to 
report that the 64 expression also agrees with that from LPT. The possible (ii, ii, 23, ii)'s 
are (0,0,0,1), (1,0,1,0), (0,2,0,0), (2,1,0,0) and (4,0,0,0). Note that the number of 
elements in {i p } n is increasing faster than n. 

The central result of this article is Eq. (|22"1) , which is a closed form expression for 
5 n , obtained for the first time, and expressible in terms of only zeroth order and on-shell 
informations. We have also explicitly checked that the results for n < 4 agree with those 
from LPT. Admittedly, the results are still not as simple as one would like to see. There 
are two sources of complexity of the results. Firstly, the cardinality of {i p } n is growing with 
n. It is easy to see that combinatorically this is the same as the number of inequivalent 
irreducible representations of the symmetric group S n [|17| . Graphically, this is equal to the 
number of ways of drawing different Young diagrams with n boxes. The other complication 
is that each f n , expressed as iterated J[ ■ ]'s, breaks down to many terms of the form J[ • ]'s. 
As a result, we expect the higher order results to grow more and more messy. 

Still, this drawback is compensated by another nice feature of the scheme, namely the 
freedom to choose the unperturbed system. Since only on-shell information is needed in 
our scheme, it is much easier to find a "unperturbed" state which approximates the exact 
system than in the case of the Rayleigh-Schrodinger scheme, in which off-shell information 
is also needed and one needs to solve the "unperturbed" system completely. As discussed in 
Ref. |T(| , one expects 8 n converges quickly as long as "Levinson criterion" (the unperturbed 



and perturbed system should have the same number of bound states) is satisfied. Hence the 
first several terms in the perturbation series should suffice to give a good approximation of 
the exact 5, and the higher order terms are seldom needed in practice. 

In conclusion, we have obtained the perturbative correction to the scattering phase shift 
of any effectively one-dimensional Schrodinger problem. The result can be expressed in 
closed form and only on-shell information is needed. With a suitable choice of the "un- 

8 



perturbed" system, the perturbation series converges rapidly. We have also demonstrated 
explicitly that the our results agree with the conventional LPT results up to the fourth order 
in the expansion series. 
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